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Abstract 

We introduce a new class of quantum models with time-dependent Hamiltonians 
of a special scaling form. By using a couple of time-dependent unitary transforma- 
tions, the time evolution of these models is expressed in terms of related systems 
with time-independent Hamiltonians. The mapping of dynamics can be performed 
in any dimension, for an arbitrary number of interacting particles and for any type of 
the scaling interaction potential. The exact solvability of a "dual" time-independent 
Hamiltonian automatically means the exact solvability of the original problem with 
model time-dependence. 
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Dynamics of quantum systems, governed by time-dependent Hamiltonians, has at- 
tracted much of attention for a long time. A relevant progress has been made mainly in 
the study of one-dimensional time-dependent harmonic oscillators^~^° which have many 
applications in various areas of physics (see e.g. Refs. 11-13). Lewis-*^ and Lewis and 
Riesenfeld^ (LR) have introduced for these systems a quantum-mechanical LR invari- 
ant, and derived a relation between the invariant eigenstates and exact solutions of the 
corresponding time-dependent Schrodinger equation. Using the LR invariant method, 
exact wavefunctions have been obtained for harmonic oscillators with timc-dcpcndcnt 
frequency,^"^ time-dependent mass and frequency,^'^ and linear driving terms.* At present, 
the exact solution can be, in principle, constructed for a general one-dimensional time- 
dependent Hamiltonian of N coupled quantum oscillators. The addition of a singular 
inverse quadratic potential does not break the exact solvability of the time-dependent 
harmonic oscillator, as has been shown by combining the LR invariant method with a 
unitary transformation in Refs. 14-17. 

In this paper, we present a new vast array of quantum models with time-dependent 
Hamiltonians of a special scaling form. By using a couple of time-dependent unitary 
transformations, time evolution of these models is expressed in terms of related systems 
with time-independent Hamiltonians. The mapping of dynamics can be performed in any 
dimension, for an arbitrary number of interacting particles and for any type of the interac- 
tion potential. The exact solvability of the time-independent Hamiltonian automatically 
means the exact solvability of the original problem with model time-dependence; other- 
wise, dynamics induced by the time-dependent Hamiltonian can be studied by simpler 
techniques known for time-independent Hamiltonians. 

We first restrict ourselves to the case of one particle situated in a d-dimensional space of 
points r = {xi,X2, ■ ■ ■ , x^)- The proposed time-dependent Hamiltonian has the following 
scaling form: 




(1) 



where p = — i/iVr is the momentum operator, V is an arbitrary one-body potential and 
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the prime denotes the differentiation with respect to the argument. The normahzation of 
parameters was chosen such that at initial time t = the Hamiltonian takes the standard 
form 

H(t = 0) = |^ + nr). (2) 

The function f{t) is assumed to be real continuous function of time in order to ensure the 
hermiticity of the Hamiltonian. To simplify the notation, without any lost of generality 
we shall assume that 

/(O) = 1, /'(O) = e (e real and nonzero) (3) 
and rewrite the Hamiltonian (|l|) as follows 
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2m fit) V/W, 
Two choices of the function f{t) are of special interest. The first one, 

/(t) = (l + 2et)V2 (4) 

with e > in order to avoid singularities in time, corresponds to the constant particle 
mass. The corresponding Hamiltonian (|l|) takes the form 

«'*) = £ + rT25^((TTl5)vj)- 

The classical version of this Hamiltonian has been introduced in Ref. 18 in connection 
with the study of adiabatic propagation of distributions, with e taken as the adiabatic 
slowness parameter. The second choice is 

/(t) = exp(et), (6) 

and the corresponding Hamiltonian reads 

Uft) = e^''^ + V(e-''r). (7) 

In one dimension and when V{x) = the model (|^) is known as the Caldirola- 

Kanai oscillator. Its exact quantum states are known, even in a more general case 
with the singular inverse square potential V{x) = 77100^x^/2 + k/x"^. 



The evolution of the system with Hamiltonian (1') is governed by the time-dependent 
Schrodinger equation 

i^^^(r,t)=H(t)^(r,t), (8) 
with the initial condition for the wavefunction 



^(r,t = 0) = ^{r). 



(9) 



In the next two paragraphs, we apply successively two time-dependent unitary transfor- 
mations, 

{H(t),^(r,t)} =^ {H,,Mr,t)} =^ {H2,^2(r,t)}, 

to go from the original time-dependent Hamiltonian H(t) to the new time-independent 
Hamiltonian H2, via an intermediate one Hi. 
We start with the unitary transformation 

1 



t' 



ln/(t), 
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J = 1,2, 
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It is easy to verify that it holds 

d_ 

dt 
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r' ■ Vr' 

J = 1,2, 



The transformation (10) converts the evolution equation to the new one 



V'i(r,t) 
Hi 



HiV^i(r,t), 

er ■ p + r(r). 

2m 



The original wavefunction in (H) is expressible in terms of ipi as follows 
and the initial condition fBI) now reads 



(10a) 
(10b) 



(11a) 
(lib) 



(12a) 
(12b) 



(13) 



7/>i(r,t = 0) = ^(r). 



(14) 



The time- independent Hamiltonian Hi in (12) can be reexpressed as follows 



Hi 



p — mer) H — e (p ■ r — r ■ pj 
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Since 



p ■ f = f ■ p — ihd, 
the evolution Eq. (12) can be written as 



Using the unitary transformation 



1 2 
— (p-mer) V^i(r,t) 



+ 



V{r) — -me'^r'^ 



V^i(r,t). 



(ie , ime 2 



V^i(r,t) = exp ( -yt+ — r 1^2(r,t) 



we finally arrive at 



d 

ih—ip2{Y,t) = H2^2(r,t), 



p 

H2 = — + 

2m 



l^(r) — -me^r^ 



The initial condition for ip2 is implied by Eqs. (|T4D and (|18D as follows 

^2(r, t = 0) = V^(r) exp ( 



(15) 



(16) 



(17) 



(19a) 
(19b) 
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The original wavefunction if) is obtained in terms of ^2 by combining relations (|T^ and 
(0), 



^(r,t) 
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' ■,-ln/(t)l. 
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Using the above mapping of dynamics, standard techniques able to deal with time- 
independent Hamiltonians (see e.g. Ref. 22) can be directly applied to time-dependent 
Hamiltonians. The mapping becomes even more appealing when the resulting time- 
independent problem with potential 



is exactly solvable, because this automatically means the exact solvability of the original 
model time-dependence. It is clear that V^(r) must be expressible as an integrable potential 
plus the harmonic-oscillator potential In one dimension, there exists an infinite 

chain of exactly solvable "reflectionless potentials". ^^'^^ In higher dimensions, there are 
standard models which admit the exact solution. We briefly discuss the simplest example 
- the isotropic d-dimensional harmonic oscillator with V{r) = mu;^r^/2. For the special 
choices (||) and (|) of the function f(t), the respective time-dependent Hamiltonians 

H(t) = p1 + Im f -^^) ' (22) 
2m 2 Vl + 2et/ ^ ^ 



and 



H(t) = e^^*^ + -e-2^WV2 (23) 
2m 2 

are transformed to the same time-independent Hamiltonian 

H2 = — + -m(uj^ - e^y. (24) 
2m 2 ^ ^ ^ ' 

The time evolution associated with this Hamiltonian can be obtained by using the method 
of the separation of variables. Without going into details, there are three regimes: iiuo > e 
one expands the wavefunction %l)2 (and, consequently, ip) in Hermite polynomials, if u; = e 
in plane waves and if u; < e in hypergeometric functions. 

The extension of the formalism to the case of an arbitrary number of particles is 
straightforward and we only write down flnal formulae. For N particles, the Hamiltonian 
(1') generalizes to 



H(t) 



(25) 



-2m, f{t) \\f{t)\, 

where /(O) = 1 and /'(O) = e, and the interaction function V [with particle coordinates 
uniformly scaled by l//(t)] involves all possible one-, two-, . . ., A^-body potentials. The 
solution of the time-dependent Schrodinger equation 

i^^V^({r,},t) = H(t)V^({r,},t) (26) 
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with the initial condition 



V'({r,},t = 0) = V'({r,}) 



reads 



V'({r,},t) 



exp 



— > rrii 



1p2 



Xm. 



ln/(t) 



(27) 



(28) 



Here, the wavefunction ■02 evolves under the action of a time-independent Hamiltonian, 



i^^V'2({r,},i) = ll2M{rjht), 



H2 = + 



N 



1 



2 2 



(29a) 
(29b) 



with the initial condition 

V'2({r,},i = 0)=V'({r,})exp 



ie ^ 



(30) 



As the simplest exactly solvable example, we mention a system of N coupled d-dimensional 
harmonic oscillators, which is solvable in the transformed picture with the time-indepen- 
dent Hamiltonian by the ordinary normal-mode technique. 

The couple of time-dependent unitary transformations presented in this paper can be 
applied directly to the quantum Liouville equation for the density matrix, extending in 
this way the treatment also to the propagation of distributions. 

In conclusion, we believe that the results presented will enable one to answer some 
problems concerning the evolution of quantum systems under a time-dependent Hamil- 
tonian, for instance, to predict a qualitative change in the analytic structure of leading 
corrections to the ideal adiabaticity when going from a few degrees of freedom to a dis- 
sipative system with infinite degrees of freedom. For classical adiabatic processes in one 
dimension, the corresponding exact analysis has been done in Ref. 17. 
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